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Abstract. Shift registers/Primitive polynomials find applications in var- 
ious branches of Mathematics, Coding Theory and Cryptography. Matrix 
analogues of primitive polynomials do exist. In this paper, an algorithmic 
approach to generating all such matrices over GF(2) has been presented. 
A technique for counting aU such n x n matrices over GF(2) is also pre- 
sented. The technique may be easily extended to other finite fields. 
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1 Introduction 

Shift registers have been widely used in the later half of the last century to 
generate sequences of O's and I's having certain other properties. A few of them 
are mentioned over here in brief. For a detailed account, one may refer to pQ. 

1. Secure and limited-access code generators 

(a) Encipherment 

(b) Privacy encoding 

(c) Multiple address coding 

2. Efficiency code generators 

(a) Error-correcting codes 

(b) Signals recoverable through noise 

3. Prescribed property generators 

(a) Prescribed period generators 

(b) Prescribed sequence generators 

4. Mathematical models 

(a) Random bit generators 

(b) Finite state machines 

(c) Markov processes 

In PP, it has been shown that each state of an n-tube shift register can 
be thought of as an n-dimensional vector. The shift register acts as a linear 
operator and changes each state into the next. This linear operator operating 
on n-dimensional vectors is represented by a n x n matrix. The characteristic 



equation of this matrix is essentiaUy the same as the characteristic polynomial 
of the shift register. In this paper, only matrices are dealt with. 

The matrices under consideration belong to the last two application areas. 

All the n X n matrices of highest order over GF(2) have been counted with 
detailed description of each step. Begin with the property of the characteristic 
polynomial of these matrices and note that these polynomials are primitive poly- 
nomials. Also note that these matrices are the primitive elements of the general 
linear group GLn{GF{2)), conjugate to each other and belong to an equivalence 
class where the the equivalence relation is conjugacy. The cardinality of this 
equivalence class is calculated from elementary abstract algebra. Once all the 
above steps are complete, it remains an easy task to calculate the total number 
of such matrices. A method for generating all such matrices is presented next. 

2 Preliminaries: A few definitions and results from 
Abstract Algebra 

At the outset, a few results are presented. For the proof and other details, one 
may refer to |ll2IHj 

Theorem 1. The maximum order of an n x n matrix A over GF{2) is 2" — 1. 

Theorem 2. An n x n matrix A over GF{2) has maximum order 2" — 1 iff its 
characteristic polynomial is a primitive polynomial of degree n. 

/O 1 

Example 1. The matrix 1 

\1 1 

has the characteristic polynomial f(x) = x'^ + x -\- 1, which is a primitive 
polynomial of degree 3 over GF{2). Hence its order is 2^^ — 1 = 7 

Theorem 3. There exist equal number of matrices corresponding to each of the 
primitive polynomials of degree n. 

Theorem 4. There are number of primitive polynomials of degree n 

where (p(n) is Euler's Totient function. 

In order to count the n x n matrices of order 2" — 1, it will suffice to count 
the matrices whose chararcteristic polynomial is a given primitive polynomial, 
say, fix). 

Since the characteristic polynomial of A is a primitive polynomial, the min- 
imal polynomial of A is also the same primitive polynomial, f{x). 
Definition: Companion Matrix 

In general, the companion matrix of a monic polynomial f{x) = ao -I- aix -I- . . . -t- 
a„_ia;"~^ -I- of positive degree n over a field is defined to be the nxn matrix 



/O ... -ao \ 
1 ... -ai 
1 ... -02 

\0 ... 1 -an-ij 

Theorem 5. All the matrices whose characteristic polynomial and the minimal 
polynomial are same and equal to f{x), (where f{x) is a primitive polynomial), 
are similar or conjugate to each other. 

The proof of the above theorem is based on the fact that all such matrices 
are similar to the companion matrix of the polynomial /(x). A direct application 
of the Rational Form Theorem proves the result easily. 

Example 2. For n = 3, there are two such polynomials. In other words, there 
are two primitive polynomials of degree 3 as follows: 

fi{x) =X^ + X + 1 

and 

/2(X) =x3+x2 + l 



'0 1 

The companion matrix of is | 1 1 

10, 



/O 1 

The companion matrix of f2{x) is 1 

\0 1 1 

The above theorem asserts that all the matrices having the characteristic (here, 
it is the same as the minimal polynomial) /i {x) are similar. Also all the matrices 
corresponding to f2{x) are imilar. 

It may be noted that the matrices we are dealing with are elements of the 
general linear group Gi„(GF(2)), the group of non-singular nxn matrices over 
GF{2). 

Definition: General Linear Group over GF{2) and its order 

The set of all n x n invcrtiblc matrices over a field F forms a group with respect 
to matrix multiplication. This group is called the general linear group of degree n 
over F and is denoted by GLn{F). The identity element of GLn{F) is the identity 
matrix and the inverse clement of A is A'^. It may be noted that GLn{F) is 
not an abelian group for n > 2. In the following treatment, we shall choose F as 
GF{2) and we shall refer to GL„{GF{2)) as G henceforth. 



The number of elements in G is given by 



n-1 

JJ(2" -2') 

i=0 

The proof of the above is easy to see: 

1. The first cohimn of a non-singular matrix must not be the vector, thus 
there are 2" — 1 many to form the first column. 

2. The i-th column must not be a linear combination of the previous i — 1- 
columns, thus there are only (2" — 2'"^) many choices. 

So 

n n— 1 

|G| = JJ(2" - 2*-i) = JJ (2" - 2') 

i=l i=0 

Let G be a group. 
Definition: Conjugacy 

If a, & G G, then b is said to be a conjugate of a in G if there exists an element 
c G G such that b = c^^ac 

Elements a,b G G related as above are called conjugate. 

Theorem 6. Conjugacy is an equivalence relation on G 

One can partition G into disjoint equivalence classes Cl{a). 

Definition: Centraliser 

If a G G, define N{a) = {a; G G\xa = ax}. It is easy to verify that N{a) is a 
subgroup of G. N{a) is usually called the centraliser of a in G. 

Theorem 7. //a G G, \N{a)\ = [g : N{a)] = 

Consider ^ as a linear operator on the finite-dimensional vector space V over 

GF(2). 

Definition: Cyclic vector 

A vector a G V is called a cyclic vector for A if the vectors A''a,k = 0,1,..., 
span V. 

Theorem 8. Let A be a linear operator on the finite- dimensional vector space 
V. Then A has a cyclic vector if and only if the characteristic and minimal 
polynomials for A are identical. 

Definition: Cyclic matrix 

If the characteristic polynomial of a matrix is the product of distinct irreducible 

factors, the matrix is said to be a cyclic matrix. 

Clearly, the characteristic polynomial and the minimal polynomial are the 
same for cyclic matrices. 



3 The counting: step by step 

Theorem 9. // the characteristic polynomial of a matrix f{x) be irreducible, it 
may be shown that any non-zero polynomial in A is non-singular. 

(The above theorem may be used since we are deahng with primitive poly- 
nomials of degree n that are necessarily irreducible) 

Proof. Let g{A) be such a polynomial. Since f{x) is irrducible, g{x) and f{x) 
are relatively prime. 

From Euclid's algorithm, 

1 = p{x)g{x) + q{x)f{x) 

for some polynomials p{x) and ^(a;). 

Replacing a; by ^ and 1 by I, the following matrix identity is obtained: 

/ = p{A)g{A) + qiA)f{A) = piA)g{A) 

since f{A) = 0, by definition of the minimal polynomial. 
Hence g{A) is invertible and therefore non-singular. 

Theorem 10. If A is a cyclic matrix and B commutes with A, then B is a 
polynomial in A 

Proof. Since ^ is a cyclic matrix, there exists a cyclic vector a such that the set 
S = V, Av, . . . , A"~^v is a basis of V. 

Consider a matrix B that commutes with A. The vector Ba can be expressed 
as a linear combination of the vectors in S since 5 is a basis. 

Ba = oo + aiAa + . . . + o„_i^"~^a = g{A)a 

where g{A) = uq + aiA + . . . + a„_i^"~"^ is a polynomial in A (of degree 
(n-1)). ■ 

As B commutes with A, it commutes with any power of A and also with 
9iA). 

Obviously, A also commutes with g{A). 

If the linear operator B is applied to any vector {A'')a in S, one gets 
B(A'=)a = {A'')Ba = iA'')g{A)a = g{A){A'')a ' 

Finally consider a matrix C = [a\Aa\ . . . \A"^'^a] 

The above result shows that the fc-th column of BC is equal to the k-th 
column of g{A)C for all A; = 1, 2, . . . , n 

which means BC = g{A).C. 

Since C is non-singular, so B = g{A). 

The final result: 



Theorem 11. The number of non-zero polynomials of degree less than n is 
2" - 1 



Proof. Consider the (n — l)-th degree polynomial cq + cix + . . . + Cn-iA"~^. 
There are n coefficients and each of them can assume only two values, zero or 
one. 

Hence there are totally 2" — 1 non-zero polynomials of degree less than n. 

Since B is a, polynomial of degree less than n in A, \N{A)\ = 2" — 1 
So the number of matrices (say,j4) with characteristic polynomial f{x) is 
equal to 

\Cl{A)\ = [G : N{A)] = = ' = H (2" " 2') 

The above is the number of matrices per primitive polynomial. 

Since there are number of primitive polynomials, so there are totally 

TT(2"-2^)^li^ tl 

f n 

number of matrices of order 2" — 1. 



4 The Method 

Given: A primitive polynomial of degree n. Let us denote it by f{x). 
Steps: 

1. Calculate the companion matrix of f{x). Call it A. 

2. Calculate the centraliser if of A in GLn{GF{2)), which is nothing but the 

collection of all the polynomials in A of degree less than n. 

\G\ 

3. Find the coset decomposition of GL„(GF(2)) with respect to H. Let = /c. 
Let the cosets be cq, ci, . . . , Cfc_i and \H\ = \co\ = \ci\ = ... = \ck-i\ = 
m(say). 

4. Choose one matrix each from these cosets. Denote the matrix chosen from 
the i-th coset Cj as Mj for i = 0, 1, . . . , A; — 1. 

5. Calculate the k number of conjugates of A using Mi, i.e. MiAMr'^ for i = 
0, l,...,fc-l. 

It is clear that two matrices belonging to the same right coset of in G 
will yield the same conjugate of A and two matrices belonging to different right 
cosets of i7 in G will yield different conjugates of A. 

Thus, if the coset decomposition of G is known, calculating the matrices sim- 
ilar (conjugate) to A is easy. All these matrices are of maximum order. 



Illustration for n = 3 



1. Notation: The following matrix 



floo OOl Ct02 \ 

aio an 012 

020 O2I 122 / 

is represented by the integer aoo2° + aoi2^ + 0022^ + aiQ2^ + aii2^ + a\2^ + 
<i2o2^ + a2i2^ + 0222^. It may be noted that there are 512 number of 3 x 3 
matrices and they will be represented by the integers 0,1,. ..,511. 
2. Coset decomposition of GLs{GF{2)) 
Note that \GL3{GF{2))\ = 168. 



/O 1 

Consider the matrix A = 1 1 

\0 1 

This is the companion matrix of the polynomial + x + 1. 
According to our notation, A is represented as 172. Construct the subgroup 
H by the centraliser of A, which consists of all the polynomials in A of degree 
less than 3. So \H\ = 2^ - 1 = 7. H = {106, 157, 247, 273, 379, 396, 486}. The 
number of cosets of H in GL3{GF{2)) is |g^='(GF(2))| ^ 168 ^ 24. Now H 
itself is one of the cosets. So there are 23 other cosets. Let us denote them 
by co,ci, . . . ,C22. 



Co = {84, 169, 253, 346, 270, 499, 423} 
ci = {85, 161, 244, 282, 335, 443, 494} 
C2 = {86, 185, 239, 474, 396, 355, 309} 
C3 = {87, 177, 230, 410, 461, 299, 380} 
C4 = {92, 233, 181, 339, 271, 442, 486} 
C5 = {93, 225, 188, 275, 334, 498, 431} 
C6 = {94, 249, 167, 467, 397, 298, 372} 
C7 = {95, 241, 174, 403, 460, 354, 317} 
C8 = {98, 281, 379, 236, 142, 501, 407} 
C9 = {102, 313, 351, 492, 394, 213, 179} 
cio = {103, 305, 342, 428, 459, 157, 250} 
cii = {106, 345, 307, 229, 143, 444, 470} 
C12 = {107, 337, 314, 165, 206, 500, 415} 
ci3 = {110, 377, 279, 485, 395, 156, 242} 
ci4 = {111, 369, 286, 421, 458, 212, 187} 
ci5 = {114, 409, 491, 254, 140, 359, 277} 
ci6 = {115, 401, 482, 190, 205, 303, 348} 
ci7 = {116, 425, 477, 382, 266, 215, 163} 
ci8 = {117, 417, 468, 318, 331, 159, 234} 
ci9 = {122, 473, 419, 247, 141, 302, 340} 
C20 = {123, 465, 426, 183, 204, 358, 285} 
C21 = {124, 489, 405, 375, 267, 158, 226} 



= {125, 481, 412, 311, 330, 214, 171} 



Now choose one coset each from the above 23 cosets and calculate the con- 
jugate of A using that matrix. Since all the matrices belonging to the same 
coset will yield the same conjugate, choosing any one of the matrices from 
each coset will give the same result. Finally we shall get a list of 24 matrices 
which are similar (conjugate) to the matrix A. The list is as follows: 
{95, 335, 187, 485, 442, 500, 102, 142, 172, 226, 106, 204, 115, 397, 157, 355, 370, 
412, 247, 431, 253, 491, 382, 478}. 

All these 24 matrices have the same characteristic/minimal polynomial, 
which is + X + 1, and all of them are of maximum order, that is 7. 

/O 1 

Next consider the matrix B = 1 

\0 1 1 

This is the companion matrix of the polynomial + x"^ + 1. 
According to our notation, B is represented as 396. Following similar lines 
as above, we construct the subgroup consisting of the polynomials of degree 
less than three and then the 23 cosets. Here, 
H = {106, 157, 247, 273, 379, 396, 486} 
and 

Co = {84, 225, 181, 459, 415, 298, 382} 
ci = {85, 233, 188, 395, 478, 354, 311} 
C2 = {86, 241, 167, 331, 285, 442, 492} 
C3 = {87, 249, 174, 267, 348, 498, 421} 
C4 = {92, 161, 253, 394, 470, 299, 375} 
C5 = {93, 169, 244, 458, 407, 355, 318} 
C6 = {94, 177, 239, 266. 340, 443, 485} 
C7 = {95, 185, 230, 330, 277, 499, 428} 
C8 = {98, 337, 307, 461, 431, 156, 254} 
C9 = {99, 345, 314, 397, 494, 212, 183} 
cio = {102, 369, 279, 205, 171, 444, 474} 
cii = {103, 377, 286, 141, 234, 500, 403} 
ci2 = {107, 281, 370, 460, 423, 213, 190} 
ci3 = {110, 305, 351, 140, 226, 445, 467} 
ci4 = {111, 313, 342, 204, 163, 501, 410} 
ci5 = {114, 465, 419, 335, 317, 158, 236} 
ci6 = {115, 473, 426, 271, 380, 214, 165} 
ci7 = {116, 481, 405, 207, 187, 302, 346} 
CIS = {117, 489, 412, 143, 250, 358, 275} 
ci9 = {122, 401, 491, 270, 372, 159, 229} 
C20 = {123, 409, 482, 334, 309, 215, 172} 
C21 = {124, 417, 477, 142, 242, 303, 339} 
C22 = {125, 425, 468, 206, 179, 359, 282} 



Choosing one matrix from each coset and then calculating the conjugate 
using that matrix finally yields a list of 24 matrices which are similar (con- 
jugate) to the matrix B. The list is as follows: 

{244, 426, 229, 171, 334, 94, 156, 354, 99, 141, 396, 114, 492, 250, 486, 190, 207, 
111, 379, 477, 415, 375, 499, 445} 

All these 24 matrices have the same characteristic/minimal polynomial, 
which is + X + 1, and all of them are of maximum order, that is 7. 

3. The matrices corresponding to the polynomial + x + 1 

{95, 335, 187, 485, 442, 500, 102, 142, 172, 226, 106, 204, 115, 397, 157, 355, 370, 
412, 247, 431, 253, 491, 382, 478}. 

4. The matrices corresponding to the polynomial + x"^ + 1 

{244, 426, 229, 171, 334, 94, 156, 354, 99, 141, 396, 114, 492, 250, 486, 190, 207, 
111, 379, 477, 415, 375, 499, 445} 

5 Conclusion 

It is clear from the above example for n = 3 and GF{2) that one can find all the 
n X n matrices of maximum order, i.e., |2" — 1|, given the coset decomposition 
of the General Linear Group of order n. These matrices are equivalent to shift 
registers, as pointed out earlier. If one knows the entire pool of matrices that 
exhibit the primitive behaviour, one can choose any one from the pool and use it 
for an application where random behaviour is sought for. If the matrix is thought 
of as a linear operator/transformation (operating on a non-zero vector of length 
n called the state), it is immediate that the resulting sequence of states will 
comprise of all possible 2" — 1 states and also this is the cycle length of the state 
diagram. However, the order of the states in each case will be different (if one 
starts with a new matrix each time), thus providing some notion of randomness. 
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